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The Hyperpolarization- activated cyclic-nucleotide gated 
channel (HCN; //, current) has been hypothesized to facil- 
itate short-term memory storage in the prefrontal cortex 
(PFC) via a cascade of signalling events involving second 
messengers including protein kinases and intracellular 
calcium (Ca) accumulation [1]. This memory mechanism 
is effected via alterations in the dynamics of the network, 
two scales above the underlying molecular level. To 
determine the factors that enable persistent activity to 
emerge, we developed a multiscale model (molecular to 
network level) of PFC that displays the persistent activity 
associated with working memory. 

The network contained 800 cells arranged in 6 neocor- 
tical layers. Cell classes included pyramidal cells, and 
fast-spiking and low-threshold-spiking interneurons. 
Each neuron contained Na, K, Ca, and HCN channels. 
Cells were interconnected probabilistically using detailed 
anatomical data from Ml, with AMPA/NMDA, and 
GABAA synapses. Intracellular Ca was admitted via 
NMDA and L-type Ca channels. Intracellular signaling 
components which contributed to persistent activity 
included: diffusible Ca, diffusible Ca buffering proteins, 
intracellular endoplasmic reticulum (ER) compartments 
for sequestering Ca, diffusible inositol triphosphate (IP3), 
ER IPS receptors (IP3R) which release Ca from ER upon 
IP3/Ca binding, sarco/endoplasmic reticulum Ca-ATP- 
ase pumps (SERCA) which pump cytosolic Ca into the 
ER, Ca/IP3 degradation reactions. Pyramidal cells had 
HCN channels regulated by protein kinases bound to 
intracellular Ca. The model was implemented in parallel 
NEURON with reaction-diffusion [2]. 18 seconds of 
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simulation time ran over 24 Intel XEON CPUs in ~7 
minutes. 

We assessed influences between intracellular signaling 
and network dynamics to isolate sources needed for 
emergence of short-term memory representations. We 
provided excitatory stimuli to a subset of pyramidal cells 
which represented the memory. Stimulus-induced depo- 
larization led to Ca influx which increased HCN conduc- 
tance and excitability of the stimulated population (firing 
increased from 0.5-2 Hz). Activation of the stimulated 
population lasted for 5-10 seconds and also increased 
beta/gamma (~15-35 Hz) oscillations. In contrast, non- 
stimulated populations were suppressed (firing rate 
reduction up to 50%) due to higher inhibition caused by 
the activated cells providing stronger drive to the inter- 
neurons. Knocking out specific pathways showed that 
persistent activity of activated cells was only viable in the 
presence of Ca influx, Ca regulation of HCN channels, 
and interactions between excitatory and inhibitory popu- 
lations, which led to network-derived inhibition activat- 
ing HCN. Lowering concentrations and binding rates of 
Ca buffering proteins reduced the network's ability to 
represent specific inputs, due to the higher levels of intra- 
cellular Ca saturating HCN channels in all pyramidal 
neurons. This predicts that Ca buffering proteins contri- 
bute to regulation of Ca at time-scales relevant for main- 
taining distinct short-term memory representations. We 
also used our model to demonstrate nonsynaptic plasti- 
city - initial stimuli cause Ca sequestration in the ER of 
activated neurons, allowing for later, selective modulation 
of their excitability via metabotropic glutamate receptor 
activation of IP3Rs, which releases the stored Ca. Our 
model demonstrates how electro-chemical interactions at 
the nanoscale may lead to network consequences that 
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produce persistent activity associated with working 
memory. 
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